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A fast, primary Coulomb blockade thermometer 
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We have measured the third derivative of the current-voltage characteristics, d 3 //dV 3 , in a two- 
dimensional array of small tunnel junctions using a lock-in amplifier. We show that this derivative is 
zero at a voltage which scales linearly with the temperature and depends only on the temperature and 
natural constants, thus providing a primary thermometer. We demonstrate a measurement method 
which extracts the zero crossing voltage directly using a feedback circuit. This method requires only 
one voltage measurement, which makes it substantially faster than the original Coulomb blockade 
thermometry method. 
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Coulomb blockade thermometry (CBT) |, § is a 
primary thermometry method which is suitable for cryo- 
genic temperatures in the range 20 mK - 30 K. It is based 
on the properties of the Coulomb blockade ||, |j| in onc- 
or two-dimensional arrays of tunnel junctions at tem- 
peratures where the charging energy Ec < k B T. Here 
Ec = e 2 /2C c ff, where C e ff is the effective capacitance 
of the tunnel junctions ||. In the CBT method, the 
first derivative of the current-voltage characteristics (IV- 
curve) is measured and from the properties of this curve 
the temperature can be extracted, using only natural con- 
stants and a calculable prefactor. The major advantage 
of the CBT method is the simple electrical measurement 
and the insensitivity to magnetic field. || 

In this letter we present an alternative approach which 
has the advantage of measurement speed. Using the same 
type of tunnel junction array we can measure the third 
derivative of the IV-curve. The third derivative has a 
zero crossing at a voltage which is (to the first order in 
Ec/k B T) proportional to the temperature. Thus, only 
one measurement point, the zero crossing, is needed to 
measure the temperature. With the original method, a 
full dl/dV vs V curve must be measured.^] 

Taking the expression for dl /dV [Q and deriving twice 
with respect to V the third derivative d 3 I/dV 3 can be 
written, to the first order in Ec/k B T 
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Here N and M are the number of tunnel junctions in 
series and in parallel (in the case of a two-dimensional 
array) respectively, i?x is the tunnelling resistance of one 



1 There is also a secondary measurement mode where only the zero- 
bias conductance is measured. This requires calibration of the 
tunnel junction resistance and capacitance. 



junction at voltages well above the Coulomb blockade and 
g{x) is defined by Pekola et al. H, and can be written 



(x/2) coth(a;/2) - 1 
2sinh 2 (a;/2) 



The function g"(x) becomes 



g"{x) = 



((a;/2) coth(a;/2) - 1) (3coth 2 (a;/2) 
2sinh 2 (a;/2) 



(2) 
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Eq. (0) is valid in the limit E c < k B T and R T > 
Rk = h/e 2 w25.8kf2. Lower temperatures and lower 
tunnelling resistances cause deviations which can be cal- 
culated theoretically (2[ fjj. In this paper we have taken 
into account the effects due to low temperature. The de- 
viations due to low resistance were not considered here, 
but are estimated to be less than 1% in the measurements 
presented here. 

We can calculate the voltage Vq at the zero-crossing of 
Eq. (Ty) numerically and the result is 



eV = ±2.U4Nk B T. 



(4) 



At low temperatures (where k B T approaches Ec) 
higher order corrections should be included (|, J|. This 
gives a correction to the zero crossing in Eq. ([j]), which 
is independent of temperature: 



eV = 2.1UNk B T - 0A65NE C . 



(5) 



The sample we have measured was a two-dimensional 
array of 256 x 256 tunnel junctions and each junction had 
an effective capacitance of 2.2 fF and a tunnelling resis- 
tance of 17kr2. The array was fabricated using standard 
shadow evaporation J9], |l(J of aluminium and in situ oxi- 
dation. The measurements were carried out by applying 
a DC voltage and an additional AC excitation (123 Hz) to 
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FIG. 1. The third derivative of the IV-curve, d 3 7/dV 3 as a function 
of the voltage, measured at three different temperatures. The volt- 
age Vq at the zero-crossing scales linearly with the temperature. 
The sample was a two-dimensional array with 256 X 256 tunnel 
junctions. 

the sample in series with a resistor i?b with a resistance 
of 20kil. The voltage over the resistor was measured 
with a Stanford SRS830 lock-in amplifier which locked to 
the third harmonic (369 Hz) of the excitation AC voltage. 
This signal SV^ is proportional to the third derivative of 
the IV-curve. The excitation voltage SV^, over the sample 
was measured with another lock-in amplifier which locked 
to the basic frequency (123 Hz) and the DC voltage was 
measured with a voltmeter. Additional low-, high- and 
band-pass filtering was used to improve the measurement. 

We used a pumped 4 He cryostat equipped with a vac- 
uum regulator which kept the bath at a constant pressure, 
and therefore constant temperature, during the measure- 
ments. 

To find the relation between the third derivative and 
SV^oj we can make a Taylor expansion of the IV-curve to 
the third order and arrive at the formula 



d 3 I 24 SVa, 



dV 3 R h SV 3 



(G) 



where 5Vsu and 5V U are voltage amplitudes (not rms 
values) . 

This method requires a relatively large excitation am- 
plitude, comparable to Vb, which introduces errors due to 
higher derivatives. Including fifth and seventh orders in 
the Taylor series we can write 



24 SVto _ d^7_ SVl^L SV * d7/ 
R h 6V 3 ~ dV 3 + l6"dV5 + 640 dV 7 ' 



(7) 



Knowing the amplitude of the excitation SVa,, the er- 
rors due to the higher order terms can be calculated and 
be corrected for. However, right at the zero crossing of 
the third derivative the higher order terms are quite small 
and therefore a relatively large SV U can be used without 
introducing large errors. 

Fig. |Tj shows a measurement of d 3 I/dV 3 at three dif- 
ferent temperatures. The shape of the curves follow the 
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FIG. 2. The temperature calculated from the zero crossing (Eq. ^) 
versus the temperature determined from the 4 He bath pressure. 



expected g"{x) behaviour, and the zero crossing follows 
Eq. |H| within a few percent. Fig. ^| shows the tempera- 
ture calculated from the zero crossing plotted against the 
temperature calculated from the 4 He vapour pressure [ [TTf 
from 1.6 K to 4.2 K. 

To test the principal advantage of this measuring 
method compared to the one used by Pekola et at we 
set up a feedback loop, illustrated in Fig. [| We used a 
DC output voltage from the lock-in amplifier, which was 
proportional to the 8V^ signal amplitude, as an error sig- 
nal to a PID regulator]^] The output of the regulator was 
used as the DC bias voltage and added to the AC excita- 
tion provided by the lock-in amplifier. The voltage over a 
resistor in series with the array was applied to the input 
of the lock-in amplifier, which was set to extract the third 
harmonic of the excitation frequency. After we adjusted 
the PID parameters to proper values, the DC voltage over 
the array stabilised at the voltage Vq, as defined by Eq. 

I 

As a demonstration that the voltage follows the tem- 
perature as expected, we took a time trace of this voltage 
while we adjusted the temperature in steps, by chang- 
ing the bath pressure, and the result is the graph in Fig. 
f|, where the voltage Vq is converted to temperature us- 
ing Eq. H The temperature steps are evident in the 
figure, and agrees well with the temperature calculated 
from the He vapour pressure, except for the lowest step. 
The disagreement at this step is probably due to higher 
derivatives and higher order corrections to Eq. |l Q At 
the beginning of the fourth step, at 600 s in Fig. ]4|, the P 
gain of the feedback was too large and the signal started 
to oscillate, but after reducing the gain (at around 680 s 
in the graph) the signal was stable again. Note that the 
relatively slow time response is not due to the thermome- 
ter or the measurement, but rather due to the time it 
takes to pump down the pressure in the 4 He bath. 

While Fig. || shows that this method is working, the 
precision is not very impressive in this first experiment, 



2 A PID regulator (Proportional, Integrating and Derivating) is a 
general feedback circuit with three adjustable parameters, which 
can be used with a wide range of applications. 
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FIG. 3. The feedback loop used to demonstrate our temperature 
measurement method. The lock-in amplifier generates a sine wave 
which is added to a DC voltage produced by the PID circuit. This 
voltage is applied to the sensor in series with a resistor, and the 
voltage over the resistor is sensed by the lock-in amplifier. The 
lock-in amplifier is set up to detect the third harmonic of the output 
frequency, to effectively measure the third derivative of the IV- 
curve. The DC output voltage is proportional to the amplitude of 
the detected signal and serves as an input error signal to the PID 
circuit. With the proper feedback parameters the voltage over the 
sensor array stabilises at Vg, as defined in Eq. 



with fluctuations up to 10%. However, there are several 
ways to improve the method. The use of a lock-in ampli- 
fier picks up the very weak third harmonic signal below 
the main excitation signal and the noise, but by notching 
out the basic frequency before the input we can increase 
the dynamic range of the amplifier, and get a cleaner sig- 
nal. The PID parameters can also be better optimised to 
the measurement, and continuously adjusted. 

Looking at Eq. [j], it is obvious that we can make an- 
other improvement by increasing M, i.e. the number of 
parallel junctions in the array. The signal amplitude in- 
creases linearly with M. This speaks in favour of using 
2D arrays with this measurement method. Note that we 
do not gain anything by decreasing N, the number of 
junctions in series, because in order to avoid higher order 
derivatives to affect the measurement we need to decrease 
the excitation amplitude SV^ by the same amount. Even 
though it was not done in this experiment, it would be 
natural to let SV W be proportional to the temperature, to 
compensate for the strong signal dependence on temper- 
ature (~ T~ 3 ). 

In conclusion, we have measured the third derivative of 
the IV-curve of a two-dimensional array of tunnel junc- 
tions. We show, theoretically and experimentally, that 
the zero crossing of this curve scales linearly with tem- 
perature, to the first order, and provides a primary tem- 
perature measurement. We also demonstrate the use of a 
feedback loop to create a fast primary thermometer. The 
feedback loop is possible because only one measurement 
is needed to get at a quantity which is proportional to 
the temperature. 

The Swedish Nanometer Laboratory was used to fabri- 
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FIG. 4. A time trace of th£ temperature calculated from the zero- 
crossing voltage Vb (Eq. ph while the temperature was varied in 
four steps. The dashed horisontal lines represent the temperature 
calculated from the vapour pressure of 4 He | Q |. The agreement 
is reasonable, except at the lowest temperature where higher order 
effects and higher derivatives affect the measurement. 
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